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We studied, both analytically and numerically, complex excitable networks, in which con- 
nections are time dependent and some of the nodes remain silent at each time step. More 
specifically, (a) there is a heterogenous distribution of connection weights and, depending 
on the current degree of order, some connections are reinforced/ weakened with strength 
<£> on short-time scales, and (b) only a fraction p of nodes are simultaneously active. The 
resulting dynamics has attractors which, for a range of $ values and p exceeding a thresh- 
old, become unstable, the instability depending critically on the value of p. We observe 
that (i) the activity describes a trajectory in which the close neighborhood of some of 
the attractors is constantly visited, (ii) the number of attractors visited increases with p, 
and (iii) the trajectory may change from regular to chaotic and vice versa as p is, even 
slightly modified. Furthermore, (iv) time series show a power-law spectra under condi- 
tions in which the attractors' space is most efficiently explored. We argue on the possible 
qualitative relevance of this phenomenology to networks in several natural contexts. 



I. INTRODUCTION AND MOTIVATION 

The concept of a network — denned as a sufficiently large set of nodes connected in pairs by 
edges — is potentially useful to help our understanding of the cooperative phenomena which are 
behind complex behavior in science and technology. Therefore, there has been a great interest for 
networks in physics during the last decade or so Most of these studies have focused on 

the case in which the edge between any two nodes is either present or not. This relatively simple 
situation allows one to investigate, in particular, wiring topology which, for example, has lead to 
the discovery of scale-free and small-world networks in natural and man-made systems. However, 
real networks exhibit a number of relevant qualities besides interesting topological structure (0; 0; 
iSS!)- In this paper, we are concerned with two features which could be essential to a network 
demeanor: 

(a) Weighted and time-dependent connections. Very generally, intensities and/or capacities vary 
notably from one edge to the other in actual networks. For instance, a main feature of trophic webs 
is the complexity of pattern flows along the food chains, the agents in social and communication 
(e.g., cell phone) networks exchange assets or information according to various rules and depend- 
ing on their partners, transport connections differ in capacity and actual number of transits and 
passengers, and effective ionic interactions constantly vary in condensed matter due to reactions 
as well as to diffusion and local rearrangements of ions and impurities. It is to be stressed that the 
connection weights in these cases often vary with time. There are variations of weights on a long- 
time scale. Their main purpose seems to be determining the nature and degree of heterogeneity the 
network needs for its intended function, say, computation, transport, cooperation, etc. In addition, 
although perhaps less investigated yet, weights may change on a short-time scale to improve actual 
functioning. To our knowledge, the best documented cases so far of such fast fluctuations do not 
belong to physics but to computational and neural networks. As a matter of fact, the human brain 
is the paradigm of a weighted network (0; [lol ) , and it is also clear-cut that high-level functions in 



the brain rely on fast synaptic changes during operation Consequently, as we have a main 

interest here on short-time weight variations, we shall in the following often use the language and 
refer to observations on neural and, eventually, computational networks. In any case, our setting 
is rather general and we believe that the main behavior described in this paper should apply to 
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networks in different contexts (see, for instance, Refs. (@;Q)). 

(b) Partial activation of nodes. In addition to the above — and perhaps also a further source of fast 
fluctuations — , one may argue that there is no need for a network to maintain all the nodes fully 
informed of the activity of all the others at all times. Relaxing such situation would both simplify 
the case and turn operation more economical. Moreover, there are some indications that certain 
nodes are more active than others, and that only a fraction of nodes is actually engaged at each time 
in some cooperative tasks. For example, this is the characteristic behavior of excitable media in 
which elements, after responding to perturbation, are refractory to further excitation 

HQ- This 

is interesting because such behavior could sometimes reveal to the observer as (relatively) fast time 
variations of connections as described above. The possibilit y of having reticent nodes is also a recent 
concern in computer science in relation with parallelism (14; 15), in mathematical-physics (jl6l ). 
and in neuroscience, where it has been associated with working memories jl?]; fl8l: fl9h . variability of 
neuron thresholds (]20l ) and silent neurons (BIB. On the other hand, there is evidence of partial 
synchronization in many different situations (|23l). In principle, this is a different phenomenon but 
one may argue that some of the observed partial synchronization processes, in which some elements 
do not attend to the others' mode, could be associated with the existence of silent and/or excitable 
units, the case of interest here. In any case, studying the effect of updating only a fraction of 
nodes will certainly shed light on the possible consequences of having partial synchronization in 
the network. 

The investigation of (a) in physics has only recently been initiated. As an example we mention 
the observed aging of nodes, e.g., in the networks of scientific publications where original papers 
stop receiving links after a finite time because review papers are then cited instead (|24l); see 
also, for instance, Refs.(0; d; @; Q; S)- However, studying the consequences of fast connection 
changes in biologically inspired models has already a two-decades history — see (I25I ) and references 
therein. For example, it has recently been shown that the susceptibility of a network to outside 
influence increases dramatically for excitable nodes 6) and , more specifically, under a competition 
of processes which tend to increase and decrease, respectively, the efficiency of synaptic connections 
at short times (j27l ). To the best of our knowledge, investigation of feature (b) is rarer (15|; 16: 128I: 



29 ). in spite of the fact that there is some — as mentioned above — specific motivation for it in 



several fields. Trying to understand the combined effect of these two features, (a) and (b), is a main 
objective here. We show that varying the fraction of nodes that are simultaneously active induces a 
variety of qualitatively different behaviors when the system is in a state of great susceptibility, but 
not in more general conditions. The susceptibility needed to observe the most interesting behavior 
is shown to occur under appropriate tuning of the connection weights with the network activity. 
It thus ensues that the effects of (a) and (b) are correlated with each other — which confirms a 
suspicion mentioned in the description of (b) above. Even more, it seems that the concurrence of 
(a) and (6) could be needed in some occasions in nature. As a first application, we describe here 
how a model exhibits unstable dynamics, which leads to itinerancy and chaotic behavior in a way 
that mimics both general expectations and some recent biological observations. 



II. DEFINITION OF A SIMPLE DYNAMIC NETWORK MODEL 

A full description of the network configuration requires both the set of node states or activities, 
s= {si} , and the set of connection weights, w = {wu S M} , where, i,j = 1, . . . ,N. From these 
we define a local field on each node due to the weighted action of the others, namely, hi (s, w) = 
X/j^i w ij s j- At each time unit, one updates the activity of n nodes, with 1 ^ n ^ N. This induces 
evolution in discrete time, t, of the configuration probability distribution according to the master 
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equation: 

P t+1 (s) = £T(s'- B )P t (s'), (1) 

s 

where the transition probability may be written as 

T(s-s') =5> n (x) [] r n { Si ^s[) J] 8 SiA . (2) 

x {i|»i=l} {i\xi=0} 

Here, x is an operational set of binary indexes — fixed to 1 at n sites chosen at each time according 
to distribution p n (x) , and fixed to zero at the other N — n sites. The choice ([2]) simply states that 
one (only) updates simultaneously the selected n nodes. The corresponding elementary rate is 



n {si ->■ s-) = a (si s-) 1 + (d s '. _ a . - l) S n>1 , (3) 



where a = a (s, (3) is a function to be determined, with (3 an inverse temperature parameter. 

The above describes parallel updating, as in cellular automata, for n = N or, macroscopically, 
p = n/N — » 1. However, the model describes sequential updating, as in kinetic magnetic models, 
for n = 1 or p — > 0. We are interested in changes with p G (0, 1) . In addition to allow for a sensible 
generalization of familiar cellular automata, this bears some practical interest, as indicated in the 
introduction. For example, assuming a neural network, p may stand for the fraction of neurons 
that are stimulated each cycle. There is no input on the other 1 — p, so that information from 
the previous state is maintained. This induces persistent activity which has been argued to be a 
basis for working memory (fl7l ; [l8l ; lish- Varying p may also be relevant to simulate the observed 
variability of the neurons' threshold (j2dl) and the possible existence of silent neurons I21I ) or dark 



neuro-matter (|22j), for instance. These are just examples of the fact that varying p has a great 
general interest to better understanding excitable media. 

The equations above may be simulated in a computer for different choices of p n and transition 
details. In order to obtain analytical results, however, we need to simplify the model somewhat. 
That is, we shall refer to the case in which the node activities are binary, Sj = ±1, the n nodes 
to be updated are chosen at random, so that one has p n (x) = (^) 5 x% — n) , and a in ([3|) 
is an arbitrary function of (only) f3s{hi which satisfies detailed balance. In spite of the latter, 
detailed balance is not fulfilled by the superposition T for n > 1, so that resulting steady states are 
generally out of equilibrium, which is known to be realistic (j3ol ). On the other hand, for simplicity, 
in order to be able to obtain some analytical results, we shall assume that the fields are such that 
hi (s, w) = h [tt (s) , £j] . Here, & = = ±1; p = 1, . . . , M} stands for M given realizations of the 
set of activities, or patterns, and tt = {tt^ (s)} , where 

^(s) = N- 1 Y t $s i , (4) 

i 

measures the overlap between the current state and pattern p. For N — > 00 and finite M, i.e., in 
the limit a = M/N — ► 0, the time equation 

< +1 (s) = pN- 1 ]T $ tanh \flhi [n (s) , £]} + (1 - p) ^ (s) (5) 

follows for any p. Actual applications concern finite values for both M and N, so that the limit 
a — > is not very intesreting in practice. This and other restrictions are not essential to the 
model, however; in fact, our simulations below concern more general situations, as pointed out 
when necessary. 
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FIG. 1 Evidence of chaos. Bifurcation diagram showing the stationary order parameter q (ir) , as defined 
in the main text, versus the synchronization parameter p for M = 20 random patterns, N = 3600 nodes, 
(3 = 20 and $ = —1/2. This behavior is characteristic of any $ ^ 1, and it follows indistinctly from the 
analytical solution and from Monte Carlo simulations. 



The model allows for different relations between the fields hi and the other network properties. 
The simplest case at hand for specific relations of such kind is Hopfield's |3ll ) which follows here 
for p — > and weights fixed according to the Hebb prescription, i.e., Wij = N Y^u^i^f- ^ ne 



symmetry Wij = Wji then assures Pt->oo (s) oc exp (PYli hi s i) ■ This corresponds to thermodynamic 
equilibrium and — using the neural-network argot — this is a case that exhibits associative memory. 
This means that, for high enough (3, the patterns are attractors of dynamics (jH), as if they 
would have been stored in the connections and recalled in the course of the system relaxation with 
time. Equilibrium is generally impeded for p > (j33l ). and the asymptotic state then strongly 



depends on dynamic details (|3Q: l34l). We checked that, in agreement with some indications (|2£ 



the Hopfield-Hebb network also exhibits associative memory for p > 0. However, no new physics 

emerges as p is varied in this case, and it is likely this occurs rather generally concerning dynamics 

for simple weighted networks. 

Our model may exhibit a complex dependence on p assu ming activity dependent weights. This 

is expected to occur in many cases, e.g., for excitable media (1121; Il3l). However, as far as we know, 

the only situation with time-dependent connections which is well documented in the literature 

concerns the brain. In this case, transmission of information and computations have repeatedly 

been reported to be correlated with activity-induced fast fluctuations of synapses, i.e., our u^,-' 
— ■ — i j 
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(jlll : l35l). For example, it has been observed that the efficacy of synaptic transmission can undergo 
short-time increasing (sometimes called facilitation) (36|; 1371 ; l38l ) or decreasing [depression) (|39l ; 
40 : l4ll ) . and that these effects depend on the activity of the presynaptic neuron. Furthermore, it 
has already been demonstrated that such processes may importantly affect a network performance 



271 : l42l : l43l : |44j). Likewise, it seems sensible to assume that similar short-time variations may 
occur in other networks — e.g., reaction-diffusion systems and the cardiac tissue ( jlH ) — associated 
with some efficacy lost after heavy work or with excitations, for instance. 

Motivated by all these facts, and also trying to maintain a well-defined reference frame, we 
shall assume that the connection weights are 



w. 



(6) 



where the second equality is introduced for simplicity. Here, Wij stands for some reference value and 



5 



0.75 




200 150 



P 



0.75 



Q. 0.5 



0.25 - 




<3> 



FIG. 2 Lyapunov surface. The Lyapunov exponent as a function of (3 for $ = 0.01 (top), and as a function of 
$ for (3 — 25 (bottom), as obtained from the saddle-point map ([5]). The black regions and curves correspond 
to a positive Lyapunov exponent, so that dynamics is then irrregular. The white regions correspond to a 
negative Lyapunov exponent associated with regular behavior. Note a small black, chaotic region for p < 1 
and low (3 in the upper graph. 



£j for a random variable. That is, we are assuming some "noise" on top of a previous preparation of 
the connections designed so that the network can perform some specific function. The background 
just described also suggests us to assume that the random variable in ([6]) is fluctuating very rapidly 
so that, on the time scale for the activity changes, it behaves as stationary with distribution given, 
for example, by 

p st (s, ej) = q5 (ej - *) + (1 - q) 5 {ej - 1) . (7) 

We shall further assume that q depends on the degree of order in the system at time t, namely, that 
q = q (irt) ■ For the sake of concreteness, our choices here will be that q (tt) = (1 + a)~ x ^ ( s ) 2 
and that Wij is given by the Hebb prescription. The result is that each node is acted on by an 
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effective field 



/if (s,w)=x; 



w 




(8) 



with = [1 — (1 — $>) q(n)]Wij. This amounts, in summary, to assume short-term variations 
which change the intensity or capacity of connections by an amount, either positive or negative, $ 
on the average. More specifically, one has a decreasing effect for any $ < 1, and enhancement for 
<& > 1, as far as <3? > 0, while $ < induces a change of sign. For the indicated choices of fields 
and reference weights, our framework reduces to the familiar Hopfield-Hebb case for 3> = 1. Note 
that it should not be difficult to implement the model for choices other than @ and (J7|). 

III. DESCRIPTION OF MAIN RESULTS 

Assuming ©, it readily ensues from (JS} for M = 1 that tt^ = F (ttoq; p, 3>) and that local 
stability requires that \dF/d^\ < 1, where 



Therefore, fixed points are independent of p for any but stability demands that p < p c with p c = 
2 {3/37T^ [(§-$) - (1 — <&) n^] - (3 + . The resulting situation for any <J> / 1 is illustrated 
in Fig. [H where one observes regular behavior, bifurcations and chaotic windows. This picture 
cannot occur for fixed weights, e.g., in the Hopfield case. In order to deepen on the possibility of 
chaos, we computed the Lyapunov surface from the analytical solution for M = 1. Two sections of 
this surface are shown in Fig. [2j This clearly reveals the existence of chaos above some degree of 
synchronization, namely, for p > p* (f3, 3>) > p c $) where the latter marks the onset of period 
doubling before irregular behavior. For example, the top graph shows that, for a small positive 
value of <I>, which corresponds to some slight depression of connections which occurs more likely 
the higher the current system order is, there is a region for large (3 (relatively small temperature, 
say T ~ 0.02 in our arbitrary units) and 1 > p > 0.8 for which dynamics may eventually become 
chaotic. In the same graph one may notice a tiny chaotic window for p ~ 1 and (3 ~ 7; this is 
the case identified previously by us ([45). The bottom graph, on the other hand, illustrates that 
chaos is typically an exception for positive values of it may only occur then for a rather large 
fraction of synchronized nodes (large p) near <I> < 0. On the contrary, for negative i.e., when 
the order tends to induce changes of sign of the connection intensities, it is more likely that the 
system will behave chaotically. It is also to be remarked that, inside the first, more exterior curve 
in each graph, there is a complex pattern of transitions from regular to irregular behavior as one 
changes, even very slightly the values of p, $ and (3. As one may imagine, this situation for very 
small M gets even more involved as M increases. Finally, it is noticeable the fact that chaotic 
switc hing among different patterns was recently demonstrated to occur also in the thermodynamic 
limit (|46l ). The next question is whether such complex behavior may have some constructive role 
in natural and man-made networks. 

Different types of behavior the system may exhibit are illustrated by the stationary Monte Carlo 
runs in Fig. [3l This involves three partially correlated patterns, as explained in the figure caption, 
and illustrates, from bottom to top: 

1. For p < p*, convergence towards one of the attractors, namely, fixed points corresponding to 
the patterns provided. This is revealed by the fact that one of the overlaps (the red one in 
this case) is constantly rather large, close to 1, while the others two are closer to zero (they 
differ from zero due to the built correlations between patterns). 



F(ir;p,<5>) = ptanhj/Jvr [l - (1 - *) vr 2 ] } + (1 - p) it. 



(9) 
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1 50 300 

Time (n) 

FIG. 3 Typical Monte Carlo runs. This shows the overlap as a function of time (in units of n MC trials), 
during the stationary regime after equilibration, for N — 1600 nodes, (3 — 50, $ = 0.004 and, from bottom 
to top, (0 = 0, 0.60, 0.87, 0.93 and 1.00, respectively. In this case, the onset of period doubling before 
irregular behavior is at p = p c ~ 0.5. This is for M = 3 correlated patterns (identified here with different 
colours). That is, we generated three patterns completely at random, and then replaced 20% of the digits 
in the second and third patterns with the same number of digits, and flipped digits, respectively, taken from 
the first pattern. 



2. Irregular behavior with positive Lyapunov exponent for a larger value of p. Notice that 
changes with time indicate that dynamics is now unstable and the system activity is visiting 
the different attractors, including the negative of some of them or antipatterns. 

3. A different type of irregular behavior in which, in addition to visiting different attractors on 
a large time scale, there are much more rapid irregular transitions between one pattern and 
its antipattern. 

4. Regular oscillation between one attractor and its negative. 

5. Rapid and ordered periodic oscillations between one pattern and its antipattern when all the 
nodes are active. 

The cases 2 and 3 are examples of instability-induced switching phenomena, namely, the sys- 
tem describes heteroclinic paths among the attractors, and remains different time intervals in the 
neighborhood of each of them, as it was previously observed in a related case (j4o1 ). 

An interesting fact concerning the nature of the phase space trajectory as p is varied is illustrated 
in Fig. HI This shows time evolution of the mean firing rate denned as 

1 - 

-i=i 

Three patterns (and their corresponding antipatterns) are involved here which consist of a string of 
Is, a string with the first 50% positions set to 1 and the rest to —1, and a string with only the first 
20% positions set to 1, respectively. In the course of this Monte Carlo experiment, we observed 
that the activity remains wandering around one of the patterns for any p < p*. The choice of 
pattern depends on the initial condition. For larger values of p within a chaotic window, as for the 
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FIG. 4 Monte Carlo runs that one may interpret as states of attention in the network, which illustrates the 
possible role of chaos. This shows phase-space trajectories of the mean firing rate for N — 1600, j3 = 167, 
$ = — |, and, from left to right, p — 0.384, 0.39 and 0.4. Here, p c = 0.38, and the system stores three 
patterns, p = 1, 2 and 3, as described in the main text. (These graphs involve a standard false-neighbor 
method (|47l ) with embedding dimension d e = 5 and time delay r = 5.) 



three cases shown in Fig. HI the system tends to visit the other patterns as well. In particular, 
the left-most case in the figure (p = 0.384) shows visits to the three patterns, and a trajectory 
which is structured, namely, there are many jumps between the pairs of more correlated patterns, 
and only a few between the most distant ones. Moreover, the number of jumps between the less 
correlated patterns tends to increase as p is further increased within the chaotic window. The 
figure shows that, for p = 0.39 and 0.40, even the antipatterns are visited; note that we have that 
£ 2 = — £ 2 . Increasing p further, e.g., for p = 0.6 in this specific experiment, the network surpasses 
equiprobability of patterns and, eventually, abandons the chaotic regime to fall into a limit cycle, 
where it periodically oscillates between a pattern and its antipattern. 




10000 20000 

Time (MCS) 



FIG. 5 Time series for the overlap 7r in the case p — 0.632, M = 1, j3 = oo (zero temperature), N = 3600, 
and <& = -0.048 (top) and -0.065 (bottom) showing chaotic transitions between the associated pattern and 
its antipattern. This series correspond to entropies S ~ 0.37 and 0.9, respectively. 

In order to deepen further on the nature of the chaotic switching, we have computed the 
normalized power s pec tra p (u) of the time series for the mean firing rate m. If one computes the 
associated entropy (|48l ). namely, S = — X^P^ log-Pw, it ensues a sharp minimum at S ~ 0.37 for 
$ = —0.048. The series corresponding to this minimum and, for comparison purposses, a different 
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FIG. 6 The power spectra corresponding to the two series in figure O i.e., for <I> = —0.065 and —0.048, 
respectively, for the upper and lower sets of data. The stright line here has negative slope 1.9. 

one for a much larger entropy are presented in figure El The power spectra for these two series is 
presented in figure [H This reveals a qualitative change of behavior, namely, that (only) the series 
describing a more efficient chaotic mechanism exhibit a power law distribution. We are presently 
analyzing in more detail this interesting phenomenon. However, we can already illustrate further 
the situation as in figure [3 This shows the distribution for the time intervals the network activity 
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FIG. 7 Distribution of the time intervals the network activity stays near each pattern during the computer 
experiment described in figura[4] 

spends wandering in the neighborhood of a particular attractor. 
IV. DISCUSSION 

We have described in this paper details concerning a model network in which connections 
are heterogeneously weighted and time-dependent, namely, correlated to the global activity. As 
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documented above, these two conditions occur in many natural networks. Furthermore, only a 
fraction p of nodes are active at each time, so that the rest maintain the previous state. This 
occurs in excitable media, for instance. 

A main conclusion is that, although the synchronization parameter p is generally irrelevant, 
varying p may greatly modify the system behavior under certain conditions. The necessary con- 
dition is a kind of susceptibility or sensitivity to external stimuli which greatly favours dynamic 
instabilities. It may be achieved in our example by appropriate tuning of two parameters, $ and 
j3. The latter is an inverse temperature which controls the stochasticity of the process. The former 
induces either enhancement ($ > 1) or lowering (<3? < 1) for positive 3>, or even change of sign 
for negative of the intensities of connections. This process is a very fast one — as compared 
with the nodes changes — , and it occurs more likely the larger the current degree of order is. The 
interesting behavior described in this paper washes out if the connection weights are fixed, even 
heterogeneously as, for instance, in a Hopfield-Hebb network, which corresponds here to $ = 1. 

Within the most interesting range for its parameters, our model exhibits heteroclinic trajectories 
which imply, in particular, a kind of dynamic association. That is, the network activity either goes 
to one attractor for p < p*, or else, for larger p, is capable of an intriguing programme of visits to 
possible attractors. The dynamic path followed during these visits may abruptly become chaotic, 
which seems the most relevant regime. Besides synchronization of a minimum of nodes, this requires 
careful tuning of p, (3 and 3>. That is, as suggested by Fig. [21 there is a complex parameter space 
which makes it difficult to predict the ensuing behavior for slight changes of parameter values. 
Note in this respect that figure [2] is for M = 3 patterns only, and that the corresponding picture 
greatly complicates as M is increased. 

The most interesting behavior of the network consists of switching among attractors. We 
observe regular switching in some occasions for p < p*, but chaos makes such process much more 
efficient. Therefore, our model confirms expectations |4^ ; [50I ; l5ll ) that the instability inherent 
to chaos facilitates moving to any pattern at any time, and that chaos and chaotic itinerancy 
may be the strategy of nature to solve some difficult problems |H3; Consistent with this, we 
have illustrated above a specific mechanism which allows for an efficient search of the attractors' 
space. More specifically, we observe a highly-structured chaotic itinerancy process in which, as 
illustrated in Fig. HJ modifying p within a chaotic window — which requires also tuning f3 and 
$ — one may control the subset of visited attractors. That is, increasing p within the relevant 
regime makes the system to visit more distant (less correlated) attractors. In this way the system 
may perform, for instance, family discrimination and classification by tuning p (1541 ) . On the 
other hand, the complexity of the parameter space for p > suggest that one could devise a 
method to control chaos in these cases. It is also suggested that one should pay attention to 
these facts when determining efficient computational strategies in artificial machines. Similar 
switching phenomena, in which the activity describes a heteroclinic path among saddle states, has 
already been incorporated in models which thus simulate experiments on animal olfactory systems 



(|5ll : l55l : l56l : 1571 : l58l: I59T ). Comparable oscillatory activity has been reported to occur in cultured 



neural networks (|19l ) and ecology models and food webs (|60j; [6l|; 162 ). This also seems to explain 



transitions between atmospheric patterns (|63l ; 164 ). and it is believed it could account for other 
natural phenomena as well (jEH). 

Finally, an important feature of the model chaotic itinerancy is illustrated in figures [7] and 
El This reveals the existence of power-law distributions within the regimes in which the network 
exhibits its most interesting behavior. This is the case for the power spectra of time series and for 
the time spent in the neighborhood of each attractor for appropriate values of p. This fact suggests 
that a critical condition which has been called for to explain some of the brain exceptional behavior 



65; 66; 671: l68l: l69l) could perhaps consists of a highly susceptible, unstable and chaotic condition 



similar to the one we have described for the model. The ocurrence of power-law behavior here is 
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consistent with the approaching to zero of associated Lyapunov exponents, sometimes referred to 
as edge of chaos. 
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